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The theory of the oscillating cantilever-driven adiabatic reversals (OSCAR) in magnetic resonance 
force microscopy (MRFM) is extended to describe the relation between an external magnetic field 
and a dipole magnetic field for an arbitrary location of the single spin. An analytical estimate 
for the OSCAR MRFM frequency shift is derived and shown to be in excellent agreement with 
numerical simulations. The dependence of the frequency shift on the position of the spin relative 
to the cantilever has characteristic maxima and minima which can be used to determine the spin 
location experimentally. 

PACS numbers: 76.60.-k, 07.55.-w 



I. INTRODUCTION 



The oscillating cantilever driven adiabatic reversals 
(OSCAR) technique in magnetic resonance force mi- 
croscopy (MRFM) introduced in ,l] has been used to 
successfully detect a single electron spin below the sur- 
face of a solid @. In the OSCAR MRFM technique the 
vibrations of the cantilever tip (CT) with an attached 
ferromagnetic particle in presence of a rf magnetic field 
cause the periodic reversals of the effective magnetic field 
acting on the single electron spin. If the conditions of 
adiabatic motion are satisfied the spin follows the effec- 
tive magnetic field. The back action of the spin on the 
CT causes a small frequency shift of the CT vibrations, 
which can be measured with high precision. 

The quasiclassical theory of OSCAR MRFM has been 
developed in This theory contains two important 
limitations. First, it assumes that the external magnetic 
field B ext at the spin is much greater than the dipole 
field Bd produced by the ferromagnetic particle. In real 
experiments, in order to increase the frequency shift 5u> c , 
one has to decrease the distance between the CT and 
the spin to values where the dipole field becomes some- 
times greater than the external field0]. Second, it was 
assumed in that the spin is located in the plane of 
the cantilever vibrations. Thus, the quasiclassical theory 
should be extended in order to describe both an arbitrary 
relation between B ext and Bd and an arbitrary location 
of the spin. This extension is presented in our paper. 

A single spin is a quantum object which must be de- 
scribed using quantum theory. The quantum theory of 
OSCAR MRFM has been developed in [4] with the same 
limitations as the quasiclassical theory. It was found, as 
may be expected, that the frequency shift 8lo c in quantum 
theory may assume only two values ±|<5w c | corresponding 
to the two directions of the spin relative to the effective 
magnetic field. The value of \8lo c \ in quantum theory 
is the same as the maximum frequency shift calculated 
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FIG. 1: MRFM setup. The equilibrium position of the spin 
and the cantilever with a spherical ferromagnetic particle, m 
is the magnetic moment of the ferromagnetic particle, jl is the 
magnetic moment of the spin, B ex t, B^ and Bo are respec- 
tively the external permanent magnetic field, the dipole field 
on the spin, and the net magnetic field. In general the vectors 
B^ and Bo do not lie in the x — z plane. 



using quasiclassical theory (where it can take any value 
between — \5w c \ and |<5w c |). Thus, to calculate the quan- 
tum frequency shift, it is reasonable to use quasiclassical 
instead of quantum theory. 



II. EQUATIONS OF MOTION 

We consider the MRFM setup shown in Fig. The 
CT oscillates in the x — z plane. The origin is placed 
at the equilibrium position of the center of the ferro- 
magnetic particle). Note that here we ignore the static 
displacement of the CT caused by the magnetic force of 
the spin. The magnetic moment of the spin, fl shown in 



2 



Fig. ^ points initially in the direction of the magnetic 
field Bq, which corresponds to the equilibrium position 
of CT. (See Eq. JI}). We assume now that the rf mag- 
netic field 2B\ is linearly polarized in the plane which is 
perpendicular to Bq. (Later we will generalize this to an 
arbitrary direction of polarization) . The dipole magnetic 
field Bd is given by: 



/io 3(m • n)n — m 
Air rf, 



(1) 



where m is the magnetic moment of the ferromagnetic 
particle pointing in the positive z-direction, r v is the 
(variable) distance between the moving CT and the sta- 
tionary spin and ft is the unit vector pointing from the 
CT to the spin. We define: 



r v = y/(x - x c ) 2 + y 2 + z 2 



y_ 



(2) 



(3) 



where x,y, z are the spin coordinates, and x c is the CT- 
coordinate (i.e. the coordinate of the center of the ferro- 
magnetic particle). At the equilibrium, the net magnetic 
field at the spin is, 



-Bo = B ext + B 



(0) 

d ' 



B 



(o) 3m^ , 

— =- zx, zy, z 
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Bext — (0, 0, B ext ), 



(4) 



(5) 



(6) 



where r — \/ x 2 + y 2 + z 2 . In the approximation which 
is linear in x c , the magnetic field B& changes by the value 



oiBf: 



B 



(i) 



'{G x , G y , G z ) x c 



(7) 



(G X7 G y ,G z ) = ^^ 7 ° ( z ( r2 ~ 5x 2 ), ~5xyz, x(r 2 - 5z 2 )) , 

(8) 

where (G x ,G y ,G z ) describes the gradient of the mag- 
netic field at the spin location at x c = 0: 



{G x , G y , G z 



(dB% dB\ dB\ 



V dx ' dx ' dx 



(9) 



(Note that the magnetic field and its gradient depend 
on the CT coordinate x c ). Next we consider the equation 



of motion for the spin magnetic moment \x in the system 
of coordinates rotating with the rf field at frequency u> 
about the magnetic field Bo. (The z axis of this new 
system points in the direction of Bq). We have: 



B, 



eff 



-7/Z x B eff , 

(Bi, 0, B - (w/7) - XcYsiGi cos cti) , 



cosai =Bq/Bq. 

(10) 

Here 0%, (i — x, y, z) are the angles between the direction 
of the magnetic field Bq and the axes x, y, z of the labo- 
ratory system of coordinates; and 7 is the gyromagnetic 
ratio of the electron spin. (7 is the absolute value of the 
gyromagnetic ratio). We ignore the transverse compo- 
nents of the dipole field Bd because they represent the 
fast oscillating terms in the rotating system of coordi- 
nate. Also we consider only the rotating component of 
the rf magnetic field. 

The equations of motion for the CT can be written, 



+ lo 2 x c = F x /m*, 



(11) 



where cu c and m* are the frequency and the effective mass 
of the CT and F x is the magnetic force acting on the 
ferromagnetic particle on CT. We consider the CT oscil- 
lations in the laboratory system of coordinates. Ignoring 
fast oscillating terms in the laboratory system, we obtain: 



F x = 



-fa Gi cos a>i 



(12) 



Next, we will use the following units: for time l/w c , 
for magnetic moment //g, for magnetic field lo c /j, for 
length the characteristic distance L$ between CT and 
the spin, for force k c Lo, where k c = m*us 2 is the effective 
CT spring constant. Using these units, we derive the 
following dimcnsionless equations of motion: 



ft = -fi X B ef f, 

x c -\- x c — F x , 

Beff = (Bi,0,A-/?fc), 

F x = -aPQ^s, ( 13 ) 
A = B Q - u, 

G = — [z(r 2 — 5x 2 ) cos a x — 5xyz cos a y 

+x(r 2 — 5z 2 ) cosa z ]. 
The parameters a and (3 are given by: 
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FIG. 2: The OSCAR MRFM frequency shift 6u c (z) at the 
central resonant surface (A = 0), for x > 0. The symbols 
show the numerical data, the lines correspond to the esti- 
mate 1171 for (a) y = (circles), (b) y = x/2 (squares) and 
(c) y = x (diamonds). Solid squares and circles indicate fre- 
quency shifts at the spin locations indicated in Fig. |3] In all 
Figures the coordinates x, y and z are in units of Lq and the 
frequencies are in units of w c . 



T 2 ' 



7/c, 



(14) 



Note that all quantities in Eq. ((13(1 are dimensionless, i.e. 
x means x/Lq, a means ////zb, Bo means jBq/lOc and so 
on. In terms of dimensional quantities the parameter (3 
is the ratio of the dipole frequency 7-B4 to the CT fre- 
quency to c , and the product a(3 is the ratio of the static 
CT displacement F x /k c to the CT-spin distance Lq. The 
derived equations are convenient for both numerical sim- 
ulations and analytical estimates. 



III. THE OSCAR MRFM FREQUENCY SHIFT 

In this section we present the analytical estimates and 
the numerical simulations for the OSCAR MRFM fre- 
quency shift. When the CT oscillates, the resonant con- 
dition ui = "f\B ext + Bd\ can be satisfied only if the spin 
is located inside the resonant slice which is defined by its 
boundaries: 



\B ext + B d {x c = ±A)\ =uh, 



(15) 



where A is the amplitude of the CT vibrations. For an 
analytical estimate, we assume that the spin is located 
at the central surface of the resonant slice. In this case 
in Eq.0A = 0. 

To obtain an analytical estimate for the OSCAR 
MRFM frequency shift we will assume an ideal adiabatic 
motion and put fi = Q in Eq. 113(1 . Let the CT begin its 
motion (at t = 0) from the right end position ar c (0) = A. 



Then the initial direction (i.e. at t — 0) of the effec- 
tive magnetic field -B e // relative to the magnetic field 
Bext + Bd and of the magnetic moment p, depends on 
the sign of Q: B e ff and fl have the same direction for 
Q < and opposite directions for Q > 0. Substituting 
the derived expression for /xj ~ — B*ifG/\B e ff\\Q\ into 
F x we obtain the following equation for x c : 



x a < 1 



a(3 2 g\g\ 



^B( + (pg Xc )* 



(16) 



We solve this equation as in |3J, using the perturbation 
theory of Bogoliubov and MitropolskyQ, and we find the 
dimensionless frequency shift (see Appendix): 



Siv c ~ - - 



ap 2 g\g\ 



{1- 



1 



tt sj B'f + {fig Ay 

B\ \ x j±jBj+WW 



2 B\ + (f3gA) 2 



Si 



.(.17) 



In typical experimental conditions we have 
B x « pGA, 

and Eq. ((17(1 transforms to the simple expression 



, 2apg 
ou) r = -- 



7T A 



(18) 



One can see that the frequency shift is determined by the 
ratio of the static CT displacement F x /k c to the ampli- 
tude of the CT vibrations A. We will also present Eq. l(T7|) 
in terms of dimensionless quantities: 



nAk r 



where 



Go = Gi cos (Xj 



(19) 



(20) 



Eqs. ((17(1 and ((19(1 represent an extension of the estimate 
derived in [3j . These equations are valid for any point on 
the central resonant surface and for any relation between 
B ext and Bd- It follows from Eq. 1 (17(1 that 5uo c is an even 
function of y and an odd function of x. 

In our computer simulations we have used the following 
parameters taken from experiments 0: 

w/2tt = 5.5 kHz, k c = 110 fj,N/m, A = 16 nm, 
B ext = 30 mT, oj/2tt = 2.96 GHz, w/7 = 106 mT, 
\G Z \ = 2 x 10 5 T/m, B l = 300 fjT, L a w 350 nm. 
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FIG. 3: Cross-sections of the resonant slice for z = —0.1 and 
z = — 1. The dashed lines show the intersection between the 
cross sections and the central resonant surface. The solid 
squares and circles indicate spin locations which correspond 
to the frequency shifts given by the same symbols in Fig. [2] 
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FIG. 4: a) The OSCAR MRFM frequency shift 5aj c (r p ) inside 
the cross-sectional area of the resonant slice for x > 0. The 
solid lines correspond to y — and the dashed lines corre- 
spond to y = x. Lines are 1 — , z — —0.1, 2—, z = —0.43, and 
3—, z = —1. r v = (x 2 + y 2 ) 1 ^ 2 . b) the cross-section of the 
resonant slice z — —0.1. The bold segments show the spin 
locations which correspond to the lines 1 on a). 



The corresponding dimensionless parameters are the 
following: 

a = 1.35 x HT 13 , = 1.07 x 10 6 , A = 4.6 x 1(T 2 , 

fix = 1.5 x10 s , B ext = 1.53 x 10 5 , w = 5.4xl0 5 . 
As initial conditions we take: 

}1(0) = (0,0,1), x c (0) = A, ± c (0) = 0. 

Below we describe the results of our computer simu- 
lations. Fig. [2 shows the frequency shift Suj c as a func- 
tion of the spin z-coordinate at the central resonant sur- 
face (A = 0). First, one can see an excellent agreement 
between the numerical data and the analytical estimate 
(|17|l . Second, as expected, the maximum magnitude of 
the frequency shift \5u> c \ can be achieved when the spin 
is located in the plane of the CT vibrations y = 0. How- 
ever, for y = x, it has almost the same magnitude \Su c \ 
(with the opposite sign of Su) c ). Moreover, for y = x 
the dependence Suj c (z) has an extremum, which can be 
used for the measurement of the spin z-coordinate. If the 
distance between the CT and the surface of the sample 
can be controlled, then the "depth" of the spin location 
below the sample surface can be determined. (In all Fig- 
ures, the coordinates x, y and z are given in units of Lo, 
and the frequency shift is in units of uj c .) 

Fig-Blshows the cross-sections of the resonant slice for 
z = —0.1 and z = — 1. The greater the distance from 
the CT, the smaller the cross-sectional area. The solid 
squares and circles in Fig.|3|show the spin locations which 
correspond to the frequency shifts indicated by the same 
symbols in Fig. [5] 

Fig. ^demonstrates the "radial" dependence of the fre- 
quency shift Su) c (r p ), where r p = (x 2 + y 2 ) 1 ^ 2 . The value 



of r p can be changed by the lateral displacement of the 
cantilever. As one may expect, the maximum value of 
\Slo c \ corresponds to the central resonant surface. The 
maximum becomes sharper as z decreases. Thus, a small 
distance between the CT and the sample surface is prefer- 
able for the measurement of the radial position of the 
spin. 

Fig, |3 shows the "azimuthal dependence" of the fre- 
quency shift Su> c (4>), where <j> — tan (y/:r) and the spin 
is located on the central resonant surface. Note that for 
the given values of z and <j>, the coordinates x and y 
of the spin are fixed if the spin is located on the cen- 
tral resonant surface. The value of <f) can be changed by 
rotating the cantilever about its axis. One can see the 
sharp extrema of the function Slo c (4>). Again, the small 
distance between the CT and the sample is preferable for 
the measurement of the "azimuthal position" of the spin. 

Finally, we consider the realistic case in which the di- 
rection of polarization of the rf field 2Bi is fixed in the 
laboratory system of coordinates. Now the angle be- 
tween the direction of polarization of 2B\ and the field 
Bq depends on the spin coordinate because the magni- 
tude and the direction of the dipole field B^ depend on 
the spin location. To describe this case we ignore the 
component of 2B\ which is parallel to Bq, and change 
B\ to B\ sin6> in all our formulas. As an example, Fig. 
demonstrates the dependence 5u c (z) for the case in which 
the rf field is polarized along the x-axis. One can see 
that in a narrow region of z the magnitude of the fre- 
quency shift sharply drops. This occurs because in this 
region the magnetic field Bq is almost parallel to the x- 
axis. Thus, the effective field Bisinf? is small: the con- 
dition of the adiabatic motion 7[i?isin6>] 2 3> \dB e ff/dt\ 
is not satisfied; and the spin does not follow the effec- 
tive magnetic field. The dashed lines in Fig. corre- 



a) 



4e~06 



2e-06 



3 

to 



-4e-06 - 



1 

y~\ 


1 1 1 1 


\ 




■ \ 




r\ \ 

\ \ 
\ * 


' / 
/ / - 
/ / 
' / 


\a 


/ / 
/ / 

/ /..•■' 
If 


\ \ 
\ 


"7 / 




/ / 

v / / 








»w / 




\ / , i 



0.2 0.4 
§1% 




FIG. 5: a) 5u} c (4>), with <f> = tan _1 (j//a;) for the central res- 
onant surface and z — —0.1 (full line); z — —0.43 (dashed 
line), z = — 1 (dotted line), b) solid line shows the cross- 
section of the resonant slice for z — —0.1. Dashed line shows 
the intersection between the plane z = —0.1 and the central 
resonant surface. The solid circle in b) shows the spin location 
4>/n — —0.1 whose corresponding frequency shift is marked 
by a solid circle on a). 



spond to the analytical estimate l|17|) with the substi- 
tution B\ — > B\ sin 9: the analytical estimate assumes 
adiabatic conditions, which are violated for small 9. 

The sharp drop of \6w c \ could be observed either by the 
change of the distance between the CT and the sample 
surface or by the change of the direction of polarization 
of the rf field. 

In any case this effect provides an independent mea- 
surement of the spin "depth" below the sample surface. 
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FIG. 6: a) 8u c (z) when the rf field B\ is parallel to the x-axis. 
The spin is located at the central resonant surface y = 0, x > 
0. Solid line are numerical data, dashed line is the analytical 
estimate 1171 . which assumes adiabatic motion of the spin 
magnetic moment pi parallel to B e ff- For a few numerical 
points indicated as solid cirlces in a) the corresponding Bo 
field is shown in b). b) solid line : intersection between the 
central resonant surface and the x — z plane. Arrows show the 
magnetic field Bo on this intersection at the points indicated 
as solid circles in a) . The absolute value of the frequency shift 
\Sui c \ drops at the spin locations where Bo is approximately 
parallel to Si (0 < 1). 
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IV. CONCLUSION 



V. APPENDIX 



We have derived the quasiclassical equations of motion 
describing the OSCAR technique in MRFM for arbitrary 
relation between the external and dipole magnetic fields 
and arbitrary location of a single spin. We have obtained 
an analytical estimate of the OSCAR MRFM frequency 
shift 5uj c which is in excellent agreement with numer- 
ical simulations. We have shown that the dependence 
Suj c on the position of spin relative to the cantilever con- 
tains characteristic maxima and minima which can be 
used to determine the position of the spin. We believe 
that moving cantilever in three dimensions, rotating it 
(or the sample) about the cantilever's axis and changing 
the direction of the polarization of the rf magnetic field, 
experimentalist eventually will enable the determination 
of the position of a single spin. We hope that our work 
will help to achieve this goal. 



Eq. (|16f) can be written in the following form: 

d 2 x c 



dr 2 



+ x c = ef(x c ), 



where r = uj c t is the dimensionless time, 

(3Qx c 



f{Xc) = 



(21) 



(22) 



and e = —a/3\Q\. 

The approximate solution of 121(1 . can be written as 





x c (t) — a(r) cos ip(r) + O(e), 



(23) 



6 



where in the first order in e, a(r) and ip{r) satisfy the 
equations: 



da , . „ , , 

— = eP 1 (a) + 0(e), 



<hp_ 
dr 



l + eQi(o) + 0(e), 



and the functions -Pi (a) and Q\{a) are given by: 



(24) 



Eq. I|28|) can be written as: 



Qx{a) = - 



:[(k 2 -l)K(k)+E(k)}, 



(30) 

where K (k) and E(k) are the complete elliptic integrals 
of the first and second kind. When k ~ 1 elliptic integrals 
can be approximated by: 



Pi(a) = — — / /(a cos i/0 sin tp dtp , (25) 
27r Jo 



i4T(fc) PS In ■ 



1 



"h, 4 



(i-fc 2 ), 



(31) 



Qi(a) = — - — / f (a cost/}) cos tp dtp- (26) 
2ira J 

On inserting the explicit expression (|22H for f (a cos ip) 
one gets: 



Pi (a) = 0, 



Qi(a) 
where 



"/ 2 (l-sin 2 y>) 
n^B 2 + (PGa) 2 Jo y/l - k 2 sin 2 tP 



{PQaf 



(27) 



dtp, 
(28) 



S 2 + {dQaf ■ 



E(k) ps 1 + - ( In 4 - - ^ 1 '- 2 
W 2 V 2 



-- a-n (32) 



In the first order approximation the frequency shift is: 



du; c ~ eQi(a) = ======= {1- 



1 



B 2 



2B? + (/30a) 



In 



V^ 2 + ggg 
Si 



(33) 



(29) In the approximation am A, one obtains Eq. l|17jl. 
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